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KINETIC THEORY OF INHOMOGENEOUS SYSTEMS 


John C. Price 

Laboratory for Theoretical Studies 
Goddard Space Flight Center-NASA 
Greenbelt, Maryland 

A theory is developed which treats the coupled equations of the 
various hierarchies as simultaneous equations in time. This 
scheme proceeds by successive approximations rather than a power 
series expansion in the small parameter ^[ n ^ d 3 ] -1 in a plasma, 
nr 0 3 in a Boltzmann gas j. The theory is suitable for the derivation 
of equations for non-uniform and force driven systems. Examples 
are given for a plasma and a Boltzmann gas. 



I. Introduction 


Although the goal of kinetic theory is the description of general non- 
equilibrium systems, most present work is directed toward obtaining small 
corrections to the behavior of infinite homogeneous systems. However for an 
inhomogeneous system the domain of validity of such corrections becomes 
smaller and smaller as their accuracy increases. Thus while the first ap- 
proximation (to the collision integral) may be presumed accurate when L » r Q 
or k d , where L is a characteristic macroscopic length and r 0 or \ a typical 
interaction length, the next correction will be quantitatively significant only 
for L > k , .. , etc. Since spatial gradients and external forces are 

mean tree path 5 sr o 

generally used to create non-equilibrium systems, the experimental verifica- 
tion of the theories will prove difficult. 

In addition the assumptions needed to obtain kinetic equations (Section III) 
are weak enough so that equations for non-uniform and force driven systems 
should be easily derivable. The fact that they cannot be derived indicates 
that some aspect of the problem has been overlooked. 

Finally it appears that most expansions in kinetic theory diverge, so that 
the mathematical simplifications must be underlaid by a physical error. 

In such a case one has reason to reexamine the basic procedures of kinetic 
theory, and to seek a first order theory (the first approximation to the collision 
term) which is as simple and general as possible. 
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Most modern work in this subject is based on the formulation of Bogoliubov , 1 
who set up the hierarchy of equations and a procedure based on expansion in a 
small parameter for solving them. In the case of a Boltzmann gas (Section VI) 
Bogoliubov's procedure breaks down for an inhomogeneous system because the 
boundary conditions (at r = o°) which he imposed are not satisfied, and are not 
relevant to the physical problem. 

His theory with some modification can be applied to yield a first order 
plasma kinetic equation for an inhomogeneous system, because three particle 
effects (shielding) are included in the plasma theory. However in this case 
the weakening of correlations described by Bogoliubov has been omitted, ^ so 
that the theory is unsuitable for non-uniform systems. 

At higher orders Bogoliubov's procedure corresponds to the calculation of 
correlation functions in successive intervals of time, which disagrees with his 
statement of the time scales involved (Section III). Nevertheless we believe 
that Bogoliubov obtained the first order theory almost correctly, and in the 
present work we attempt to correct the difficulties which appear in the higher 
orders. 

Sections II and IH, which are non mathematical, discuss the motivation for 
the work, the significance of the small parameter, and the role of time. The 
statement of the procedure to be used appears in Section IH-D. 

In Section IV the first order plasma theory is developed, and in Section V 
the second order theory is compared to the quasilinear theory^ and the recent 
work of Dupree.^ 
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Section VI contains a discussion of a Boltzmann gas, while Section VII con- 
tains a brief application to an equilibrium plasma. 

In the conclusion the general outlook for the theory is discussed. 

II. The Purpose of Kinetic Theory 

In a mathematical sense the description of a complex physical system 
(N ^10 23 particles) by one or several continuous functions obeying relatively 
simple equations is not justified. The equations of motion are known (ignoring 
quantum mechanical effects, and in the present treatment, radiation) so that in 
principle there is no need for further approximations. 

In fact, of course, one cannot solve these equations for the motion of the N 
particles comprising the system. In addition the necessary boundary condition 
(typically the initial state of the system) cannot be given by experiment, so that 
even a formal solution to the equations of motion is not useful. 

Thus our theory is necessarily statistical. We believe that the justification 
for our procedures will ultimately come from statistical mechanics, but this 
justification is lacking at present, except in the case of thermal equilibrium. 

It follows that the theory is completely ad hoc until experimental evidence is 
available. In this respect verification of the Boltzmann and Vlasov equations 
represents the major evidence that the theory is headed in the right direction. 
For the same reason we attempt to drive theory toward experiment by con- 
sidering non-uniform and force driven systems, as they typically represent non- 
equilibrium experimental situations better than homogeneous field free systems. 
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One may conceive of experimental verification of the theory at many levels. 

A small detector /AV ~ \ 3 , .. , AV ~ \ n 3 , ) may measure the directed 

flow of particles in order to measure the one particle distribution as a function 
of time . A pair of such detectors spaced closely together may be used to ob- 
tain the evolution of the pair correlation function, etc. 

In fact such detailed measurements are virtually impossible, so that only 
crude comparisons with the predictions of kinetic theory may be achieved. For 
this reason we treat basic physical principles as internal constraints on the 
theory, even though direct experimental verification may prove difficult. In 
the present work the conservation laws of particle number, momentum and 
energy will be used as a test of the theory. 

III. Fundamental Assumptions 

A. Mathematical Framework 

Since the procedures to be used in setting up a kinetic theory are a matter 
of choice we must be careful to justify the equality sign in any equations we use . 
For this reason we work with a hierarchy of equations which can be derived 
rigorously. The Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy is natural 
for a Boltzmann gas (Section VI), while the Klimontovich-Dupree ^ equations 
are more convenient for describing a plasma (Section IV). We also work with 

O 

another (quasilinear) hierarchy (Section V) first mentioned by Dupree, because 
it is mathematically simple. We believe that these latter equations are unsuitable 
for describing real systems because of the absence of source terms in f , but they 
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permit a direct comparison with quasilinear theory and its extension by Dupree. ^ 
It should be emphasized that the functions satisfying the above equations are not 
the quantities of physical interest. Each closed set of equations obtained by 
setting the (n + l)th correlation function equal to zero is fully time reversible, 
while we seek equations describing an irreversible approach to a limiting state, 
e.g. thermal equilibrium. 

In what follows we justify the continued use of the equality sign by keeping 
our equations formally identical to those of the respective hierarchy. We ob- 
tain successive approximations to a set of correlation functions, where the dif- 
ference between the exact functions and the approximate functions remains in 
the equation and may, in principle, be determined. 

B. The Use of the Small Parameter 

We consider each of the three hierarchies as a time reversible (and hence 
improper) approximation to a different set of equations (a kinetic hierarchy) 
which describes the quantities of physical interest. The arguments that follow 
apply to the physical quantities rather than the solution by direct integration of 
the various time reversible equations. 

We may begin by specifying the relative importance of terms in each 
equation by scaling velocities, lengths, times, etc. to those which are believed 
to dominate the physical behavior. We do not carry out the procedure here, 
but simply state that it leads to the appearance of a small parameter. This 
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parameter (to be called fS) is essentially nr 0 3 in a Boltzmann gas, and (n^ Debye 3) 1 
in a plasma. In most physical situations the numerical value of p is 10 ~ 2 or less, 
but it may be somewhat larger without destroying the significance of the theory. 
We choose to work with unsealed quantities, and insert a coefficient /? = 1 where 
the small parameter appears in the scaled equations. The significance of p in 
the theory we wish to set up is the following: 

We demand that the sequence of functions obtained by including more and 
more physical effects (i.e., additional terms in p in the equations) be convergent. 
Thus 


T(/3 n + 1 ) - ?(/?") 

V(/3") 


0(/3 n ) 


This severe requirement is weakened by several conditions which cannot 
be avoided. 

1. The requirement cannot be satisfied over all of phase space. For ex- 
ample the ordering in /3 breaks down for short distances | r. - r . | 

- q 2 /mv tbermal in a plasma, and for large distances I r. - r. | ~^ mean freepa th 
in both a plasma and a Boltzmann gas. Thus convergence is strictly 
required only in the region specified in the original ordering. 

2. For short times (t < t collision ) the ordering procedure is meaningless. 
Here we use procedures discussed in Part C. 

3. It is possible that solutions may diverge after long times. This is un- 
important provided the system effectively reaches an end state before 
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the divergence occurs. We find asymptotic behavior by special pro- 
cedures such as transport theory (e.g., Chapman- Enskog^ theory) rather 
than by direct integration in time . 

Although the demand for convergence does not specify a mathematical pro- 
cedure, we believe it eliminates the one most commonly used. In general a 
perturbation expansion in fi maps the flow and scattering of particles onto 
higher and higher order terms, while in fact they are moved about within the 
same function. The expansion breaks down after the displacement of a sub- 
stantial number of particles in some region of phase space. Because expansions 
may be expected to diverge, we do not expand . Instead, we truncate at a given 
level consistent with the "small" statistical effect in the next higher equation, 
and attempt to solve the resulting equations exactly. At a given level a perturba- 
tion expansion may be convergent for numerical estimates, but the expansion 
must be reappraised at each higher level of approximation. 

Obviously convergence does not determine a mathematical procedure, it 
only puts a "boundary condition" on the techniques that may be used. We believe 
that the physical behavior is represented by an intricate mathematical de- 
pendence on the small parameter, so that a statement of a general mathematical 
procedure is not possible. 

C. The Significance of Time 

Our conclusions regarding time are similar to those of Bogoliubov. 1 Our 
statements apply to quantities which we might conceivably measure in the 
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laboratory, for which we wish to obtain equations. If we consider the behavior 
of a physical system which is set in motion at some instant of time, then the 
evolution may generally be broken into three phases. 

1. In a time of the order of a collision time t c an arbitrary initial state 
will relax so that the correlation functions become substantially func- 
tions of the one particle distribution f, and hence of one other . Some 
effects of the detailed initial conditions may persist (e.g., in the case 
of plasma instabilities); these must be included in the theory. During 
the initial relaxation no general procedure (including ordering) is valid, 
and kinetic theory is entirely an initial value problem. 

In the mathematical theory integrals over the product of a potential 
and a correlation function appear frequently. If the initial conditions 
are reasonable then these integrals may reach their asymptotic form 
muchmore rapidly (t - r 0 /v thermal or t - \, ebye /v thermal ). This is 
the relaxation time described by Bogoliubov. 

2. Following the initial relaxation f decays roughly as <9f/dt = l/t c (f - f 0 ), 
i.e., rapidly at first and then more and more slowly while approaching a 
local end state f 0 . In the early stages the separation of times 1 and 2 
breaks down and the theory we develop is less valid. This is not im- 
portant, as we are generally interested in the asymptotic behavior of f , 
rather than the precise way it gets there. The aim of the present paper 
is to find equations valid during this second regime. These equations 
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should describe the final relaxation of the correlation functions as they 
become functionals of f, and the subsequent evolution of f. 

3. The third time is often given by L/v thermal » but it is better described as 
the decay time for a macroscopic state, e.g., a thermal gradient or a 
gaseous shock. In this regime a description by transport theory should 
be adequate. 

The second and third time intervals are described by the solution of the 
kinetic equation, while only the first time interval is relevant for the derivation 
of kinetic equations. In particular we consider the end of the first time interval, 
as the system approaches asymptotically the kinetic regime. Difficulties (e.g., 
the breakdown of ordering) early in the first phase are avoided simply by con- 
sidering t always in this asymptotic limit. 

In a plasma some effects of initial conditions may persist in this limit. If 
a general equation might be written 

(d \ <9W n 

+ H ) v + *<* = 0) = £(t)i>(t = 0) 

then the term on the right is permissible, providing the operator £( t ) eventually 
carries £(t ) $(t = 0) to 0. The term 3>(t =0) (9*P/<9x is improper as it continues 
to affect the evolution indefinitely. This is contrary to our demand (relaxed 
perforce in Section VI) that the evolution of the system depend only on f after 
sufficient time . In addition this term must vanish if the evolution is to be 
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independent of the origin of time. This reasonable requirement also implies the 
condition £(t ) £(r) = £( t + t). 

In what follows all equations are treated as simultaneous equations in time. 

D. Formal Procedure 

The procedure we advocate is one of successive approximations. For 
simplicity we discuss the first order statistical correction to the fi = 0 equations. 

1. The equations of the hierarchy are integrated for the case fi = 0. It is 
important to keep the streaming (homogeneous) solutions for the cor- 
relation functions, although these terms frequently may be neglected in 
the solution of kinetic equations. The resultant expressions for the 
correlation functions V(/3 ~ 0)are then substituted into the small (/3) 
terms in the hierarchy. The difference between these expressions and 
the (unknown) exact solutions remains in the equation, but is now higher 
order 'P-'P(/3 = 0) = 0(/?P), and would be considered in the second order 
theory. 

2. The small term /31 / n + 1 will, in general, contain an integral over the 
earlier behavior of lower order correlation functions 

/^ n + 1 ~ f dT k(T)' Vi(T) ••• f(r)l 

Jo 

The behavior of the lower order functions is now approximated by the 
(3 = 0 solutions such that all correlation functions are evaluated at the 
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same time t. Thus 


'M'O = £ n ( T > + { }dr 

A 

where the operator £ n is the streaming operator that carries v H n from 
time t to time r. 

3. The <P n (t ), *P n _ 1 (t ) • • • are now considered to be the exact correlation 
functions (rather than the /3 - 0 correlation functions) so that the ap- 
proximation is thrown over from the functions to the operators C. As 
we are deriving equations rather than solving them, we are approximating 
(scattering) operators rather than functions. 

4. The scattering operators involving time integrals are now evaluated. 
There is now a difficulty which appears in the first order theory, and 
presumably in the higher orders. The operators just defined are not 
mathematically proper, e.g., they may diverge. This must be corrected 
by making each operator consistent with the small (order /3) operators 
in the next higher equation. The reason for this difficulty is found from 
examination of the resulting equations: we are using the approximation 

= 0 in order to get a grasp on the problem, but the desired result is 
not analytic at ft = 0. In general this consistency is the most difficult 
part of the theory to obtain, but since it represents a small (order f5) 
effect on an operator which is already small, the theory is not sensitive 
to the exact method chosen. 
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The resultant equations represent the kinetic hierarchy including 
the first (order fi) statistical effects. 

Although we can easily generate a formal statement and notation for going 
to arbitrary order in /3, the added weight and complexity add nothing to the 
theory. Instead we describe the general procedure only so far as we are able 
to carry it out explicitly, and to the level of the first equation we cannot solve. 

We state two features that are not immediately apparent. Firstly, we are 
obtaining successive approximations for a correlation function ¥ as a function 
of lower order functions; ¥ = ¥ (¥ , , ¥ 0 , * * * , f). It follows that the 

kinetic hierarchy closes naturally at each level n. Secondly, the analytic solu- 
tion at the level ^ is required for the construction of the /3^ +1 equations. This 
illustrates forcibly that it is not possible to obtain corrections to equations 
which are already insoluble. 

Despite the fact that the theory is convergent by definition (if this cannot 
be arranged a new ordering is necessary), so that one may in principle go to 
any order, we have no proof that this is a correct procedure. The question of 
validity lies outside the theory. 

IV. First Order Plasma Theory 

Our techniques follow the methods developed by Dupree 7 in his brilliant paper 
of 1963. The equations are those developed by Klimontovich and Dupree; they 
are equivalent to those of the BBGKY hierarchy. The one particle distribution 
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for species p, f (rj, Vj, t) is normalized so that 



N m ( r i- 0 


n 


where N. t is the local particle density and is the system average density, 
f satisfies the equation 



= - 0 sr 37; • < 8t , SE > <i) 

where - (r . , v . ) , fi is the ordering parameter, and the second notation for 
the collision term is often convenient. 

The electric and magnetic fields (EandB) are determined from Maxwell’s 
equations, and will be regarded as known. We shall treat the forces as con- 
stant in time (an adiabatic hypothesis), which is valid if the correlation functions 
reach their asymptotic values rapidly compared to the variation of the fields. 10 
In our examples we consider forces which are constant in space, although they 
may in general vary over distances much greater than a Debye length. 
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The zero order (/? -- 0) solution to (1) is f (X x , t ) = e V(1) t f ^ (X t , t - o) 
where the operator e -Vt runs the particles backward on their orbits 


r(t) 


r (t), 


v(t) 


m 


y(t) 

E(r) + — x B(r) 


from the initial point r(0) = r ; v(0) = v. If the solution to the zero order 
(Vlasov-Maxwell) equations is not known it is not possible to go to first order. 

As implied by the right side of Eq. (1), we use the Coulomb approximation in 
treating fluctuations in the system. Electromagnetic effects may easily be in- 
cluded in the theory, 11 but they complicate calculation considerably. For 
present purposes any fluctuating field E is related to the fluctuation density f by 
Poisson's equation 


E(r, t) 


v * 


dX' V 


r - r 


t) 


In order to permit comparison with Section V, we write the equation for the 
general correlation function (ignoring particle species 1 ^) h n (X t , X 2 , * * • , X n ; t ). 
We define the Vlasov operator 


V(i ) 


v . 




E + 


v 

1 

c 


X 


B 


d 

dv . 
1 
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the fluctuation operator 


T(i) 


V(i) 


q_ df 

m dv . 
1 


f 


dX.' V. 


nq 


and the interaction operator 


0(i> j) 


~q 

m 


A | 


dX. V. 
] 1 


nq 


- r . 


Then h n satisfies the equation 


dt 




’) ^n-m+l ( " ’ ^n + l) + ^ 1 n+l( 


= 0 

( 2 ) 


In the sum on m the remaining coordinates are distributed in all ways such 
that each appears once, and rearrangements within a function are not distinct. 
The correlation functions h are related to those of the BBGKY (g , with 

n 11 

g 2 = "g"> g 3 = " h "> etc -) b y 




+ • • • + f W • X 2 ) A ( X , ’ X 3 ) ' • ' A ( X . ' X n) W 
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where 1 3 


Aft.X,) 




The first correction in f3 to the equation for f may be found by calculating 
h 2 to zero order. We set /3 = 0 and n = 2 in Eq. (2), and use the more 
familiar notation h 2 = (Sf8f). 


d_ 

dt 


T(l) + T(2) 


) <SfSf> 


(4) 


As pointed out by Dupree a product solution is possible. We define a 
propagation operator P(X, t) by 



P 


0 


(5) 


with the initial condition P(X, 0) -I. It is convenient to define an auxiliary 
(electric field) operator 


P(r, E, t) 


f 


-nq dX' V 



P(r\ v', t) 


(6) 


so that 5 may be written 



We first integrate along the orbits as described previously. 


P - e vt - I dr e" v < t-T) P(r, E, r) • f(r, v, r) 

This expresses P in terms of f at other positions and earlier times. However 
from the /3 = 0 approximation for f we have f(r) = e V(t-T) f(t). Thus 


P 



dr e _v ^" T) P(r, E, r) • 


d 

<?v 


f(X , 


t) 


(7) 


Operator Eqs. (6) and (7) may be solved by a Fourier transform in space 
(k) and a Laplace transform in time («), provided we use an adiabatic hypothesis 
to ignore the space and time variation of f in the transforms. Then the trans- 
forms of the operators are given by 


P(k, v, co) 


- f 

•'n 


dre ik-[r ( -r)-r] e -vr e i + SL'A 

K 


1 n .co 

dre i “ T e ik ’[ r <- T >‘ r ] 

o 



d dv(r) 
dr + d\ 



v, t) 


47t nq 

e(k, o)) 



*’ [ r '^ r )- r '] e 


i COT e “VT 


( 8 ) 


P(k, E, Cl ) ) 


jjk 

k 2 


47t nq 
e(k, w) 



dre iBT 


e ik [ r <- T H 


e 


-Vt 


(9) 
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where the generalized dielectric function is given by 


e(k, co, r, t) 




Jdv 


[ 

*0 


dr e 1 


e ik-[r(-r)-r] 


dr(r) 

<9v 


d (9v(t) 
dr + <9 v 



f(r, v, t) 


(10) 


To use the P operators we Fourier transform the operand, apply the oper- 
ator, and invert the Fourier and Laplace transforms. For example the solution 
to Eq. (4) is given by 


(8f8f | X lf X 2 , t) = p(X lt t) P(X 2 , t) (SfSf | X 2 , t = 0) (11) 

This solution is improper in several respects. For example one part of the 
solution is given by (SfSf | X t (-t), X 2 (-t), t = 0), which does not vanish as 
| r. -r. | -<». This defect is caused by the omission of three particle effects, 
which tend to smear out the orbits and prevent the particles from being cor- 
related to | r . - r . | -oo. This represents non-analytic behavior in /3, and should 
be corrected in the second order theory. 

A second defect arises in the boundary conditions which should be applied 
to the correlation functions. The difficulties of using the Fourier transform 
illustrate this, for in general (§f8f) is not simply a function of r . - r . . One 
cannot seek the asymptotic behavior (in time) of wave effects^ for they cor- 
respond to the propagation over large distance of the initial fluctuations, which 
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may be determined only by specifying a boundary condition. Once again we be- 
lieve that the proper resolution should be found in the treatment of three particle 
correlations. 

Despite these difficulties we may find a first approximation to the kinetic 
equation, for the potential q/ 1 r x - r 2 1 in the collision integral cuts off most of 
the dependence on | r x - r 2 | > ^ Debye . We ignore the "slow" dependence on space 
and approximate 

(SfSf | kj, k 2 , v x , v 2 » t = 0) = f (X r t = 0) ( ' 2 ^ - S(v x -v 2 ) 8(k x + k 2 ) 


+ g 2 (kj, k 2 , Tj, v 2 , t - 0) (2 7T) 3 S(k x +k 2 ) 

The fimction f (X x , t = o) is written e V(1)t f (X x , t) in the collision integral, where 
the operator e V(1)t must be inserted in the time integrals defining the trans- 
formed (oo) P operators. If we shorten notation by defining 


r V x (k, v, 


P v (k, v, co) 
P o ( k > v * w ) 

[^l (k, V, Cc>) J 


r 


dr(r)~ 

<9v 


<?v(r) 

k ' -57- 


dr e 1 0)7 e 


,ik-[r(-T)-r] 


e - yT ^ 


( 12 ) 
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then the kinetic equation is given by 




477nqi 
k 2 e (k, ojJ 





dv 2 P 0 



X 







P v (k, 






47rnq 

e(-k, o) 2 ) 


X 



The Laplace transforms are evaluated by considering t in the asymptotic 
limit. We illustrate the procedure with a simple example. 
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A. Inhomogeneous System 


For a system with no fields the velocity visa constant and r( t ) = r + vt. 
The operators P 0 , P r andP v are given by 


P 0 - (-ioj+ik'v) 1 P f - -(aj-k'v) 2 k P v 

and the dielectric function is 


(-ico + ik *v) 1 k 
(14) 


e(k, co, r, t) - 1 + 


2 . , , df 

* p f dv k • av v r uv k - ar 

k 2 J k • V + k 2 J 


dv k 


df 


k 2 J (co - k • v) : 


(15) 


As the methods of evaluation have been discussed elsewhere 15 we simply 


give the result 


df 


dt + v i ’ dv. 


f = 


2q 2 _d_ 
77Tn dv 


; I 


dk 


kk 

k 4 


r 


J 


nq 2 dv.< 


7T8 (k • Vj - k • v 2 ) 


e(k, k • v 2 )| 


®(k, k • v 2 ) 


m ~ al;] f ( r i- v i) f ( r i- v 2 ) 


(k • v t -k • v 2 -i p) 2 


® (k, k • Vj) ! 


e(k, k • v 1 ) | 2 m f ( P l’ Vl ) 


a ®(k,k-v 2 ) j a 

^7 f ( r l’ v 2 )- f ( r l- y 2 )^ f ( r l’ V l) 


+ j( k - r i> v i) 


S (k ' v i _n k) 


A 


f ( r l> v i) + 2 [ dr' E(k ’ t ’ ^ f(r r v v t') S[k- Vl -n k (t')] 1+-^- j* 
de , J 2 I J 

d co v k) 


£w 


dv nq I dv 


i 


¥^4*', -.) 7 

dco K) 1 


(k. * 2 > V 3’ 1 “ °) 


n k -k • v 2 -i 


i |y k |) (n k -k • v 3 + i |y k |) 


(16) 
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t 


where p is a real infinite smal.o^ - + 7 k , and we have defined 


j( k < r i- v i) 


\a) k " k • v I 


277 S (k • v - O k ^ 


1 di ( r v V l) 


1 of( r i, Vj) 


m 


(k • Vj - k • v 2 - ip) 


2 m dr, 


© (k , k • 

v) 

= ±1 

Im e(k, k • v) < 

E(k, 

t\ 

t ) = exp 

2 J Tu (r) dr 


The collision integral conserves number density and momentum, but it does 
not yield a correct energy law, for the energy transport cannot be written as 
the divergence of a flux. Although the difficulty is easily traced to the spatial 
adiabatic hypothesis in the Fourier transform, we believe that a more correct 
mathematical treatment is not justified, as this will describe the flow over 
large distances of the microscopic fluctuations originally present in the system. 
Instead a treatment of three particle correlations should cut off the dependence 
of the collision term on distant processes, as well as conserve energy. 


B. Plasma in a Magnetic Field 

For simplicity we assume a plasma in a uniform magnetic field Be z al- 
though slow variation of B in space, and weak electric fields may be included 
easily We express velocities and wavevectors in cylindrical coordinates, 
taking the x direction as origin for angular variables; v fv it 6 , v z ) , 
k (k it a, k J. For physical applications it is appropriate to keep spatial 


dependence in cartesian coordinates, so our notation is sometimes mixed. The 
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particle orbit is given by 


r(t) = r + e z v z t + e x |sin(<9 + oj c t) - sin^ + e y ^ [cos0 - cos (d + w c tjj 
where 


OJ 


qB 


c me 


The P 0 , P r and P v operators are given by 


r ^ 
P 


r 


<*•> ‘ L 


e~ ik l v l/ CJ c lsin (^" a ) e in( 0-a) 




U„(a) 


_ .. «J„ 


r 


< k z Jn ( a ) + T"Jn ( a ) e X + ik l e 6 


> 


ik z J n (a) L 2 k x L, - J n (a) 


co 


kx B 




where 


n <?J n 

Lj (n, a, 0, a) = ^ J n (a) cos (6 - a) - i ~s~ sin (6 - a) 

n ^Jn 

L 2 (n, a, 8, a) =-^-J n (a)sin(0-a)+i cos (0 - a) (18) 
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The Laplace inversion and analytic continuation are now straightforward 
but the result is so intricate the evaluation should be performed after approxi- 
mations suitable to the problem at hand. A single difficulty arises, which is 
the appearance of terms oscillating at multiples of the cylotron frequency. 

These should be eliminated by time averaging over one period, for if the adia- 
batic hypothesis is not satisfied (in effect making df/d6 small) then our analysis 
is not relevant to the rapid time variation of the system. 

For purposes of approximating the k integral one must know the dependence 
of the dielectric function a. However this lies within the linear theory and will 
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not be discussed here. Even in the case of general fields leading to particle 
drifts the collision integral conserves particle number and canonical momentum. 
The transport of electrostatic energy is not described properly so that a more 
thorough study of wave effects is needed. 


V. Toward Second Order Plasma Theory 

Our treatment in this section is very incomplete because of the great com- 
plexity of the mathematics. We simplify the problem by working with the quasi- 
linear hierarchy described by Dupree. 8 The equations are formally identical to 
those of the preceding section, but the correlation functions (w n ) have no delta 
function between spaces; they are non-singular. One obtains the same equations 
by dropping from the BBGKY equations all terms in V0, while keeping integrals 
over such terms. Because the w n equations have no sources in f the correlation 
functions do not relax to become functionals of f . In this respect the system 
never becomes "kinetic" for knowledge of f alone is not sufficient to determine 
the subsequent evolution. 

Using the operators defined in Section IV, we have 


d 

df + 


n 

L 

i= 1 


T(D 


w 


J2 ° (i> n+ 

x = 1 


n- 1 

^ W ™("‘’ X i’ '") W n-m + l I”*’ X n+l) 


m - 2 


+ 


• • • X n+l)f( 20 ) 
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The streaming solution (including shielding) is given by 


W n , 0 - P n ( X l- ' ’ ' X n= t) W n “ 0 > (2 ^ 

where we have abbreviated P n (Xj, X 2 , • * • X n ; t) = p(Xj, t) p(X 2 , t) • • *P(X n , t). 
Equation 20 has been integrated by Dupree for the special case Vf = 0, i.e. for a 
homogeneous system with no forces, or else df/<9<9 = 0 in the case of a uniform 
magnetic field. We shall use his solution, despite the fact that it is not strictly 
correct in general. 1 ? 

For ft - o we have 


w 


w 


r 


dr P ( t - t ) 


n n - 1 

2Z 0(i ’ nti > TL 


w 




o - 2 


w 


“ m + 1 


(•••■ 


1 ’ 


r) (22) 


Since the source terms in the integral are of lower index an iterative solu- 
tion for w n is now possible. Instead we use successive approximations in the 
source terms, starting with the streaming solution. We approximate the be- 
havior of the w m by the streaming solution, after which the w m are treated as 
the exact functions. We then use the commutative properties of the P operators 
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acting on different spaces to write 


i t n n~ 1 

drP n(t-r) ^ 0(i. n+ i) 

i = 1 m - 2 


-r) 


P 


- l 

n~ m + 1 


(t-r)w m (t) w n _ m+1 


(t) 


P n (t>W n (t = 0)- 



t 


n~ 1 

-r)0(i, n + l)P _1 (X., t-r)F rl (X n+1 , t-r) 

m= 2 


'>• t ) W n- m +l(‘" X n+l’ f ) 

(23) 


Substitution of w ntl l into the equation for w n yields 


n n I n ~ 1^ 

^" + ZZ T(i) W " = ~ 7"! Q(i ’ n+ ] J"! (• • • X i ’ * * *) W n~ m + 1 (‘ * * ’ X n+l) 

i=l i = 1 1 m ” 2 


+ 


W 


+ 1.1 



The first w n + 1 j should be considered order 1, for while it is quantitatively 
small its effect on w n cannot be found by expansion. The terms with factor /3 
should then be evaluated from the new /3 = 0 equations. 

Because w 2 0 is the /3 = 0 solution for the source terms of w 3 , it follows 
that w 3 j is the correct /3 = 0 solution for w 3 . Thus /3(w 3 - w 3 j) is formally 
order /3 2 . For n > 3 the source terms must be iterated n - 2 times before 
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/3 (w n - w n n _ 2 ) is order / 3 2 . We pursue the general case no farther, but consider 
n = 2, and drop the /3 2 terms. 


<9t 


+ 


T( 1) + T( 2) 


) 


w- 


2 


L 

j = i 


0(j, 3) 



t - t) 0( i > 4) P 1 (X. , t - r) 


2 

x P -1 (X 4 , t-r) w 2 (-,X.;t)w 2 (-, X 4 ;t)- 0( j , 3) w 3 0 (25) 

j = 1 

Although w 3 0 depends explicitly on the initial value of w 3 , the evolution of 
w 2 should not depend on a particular choice of the origin of time. To ensure 
this we must give a consistent treatment of the initial value problem. 

Because of the asymptotic analysis (Section III-C) and equation (23) we have 


1 t-“ 3 

dT J> t-r)0(i, 4)P" 1 (X., t-r) P’ 1 ^, t-r) 


W 3, 1 (O " P 3 (t ) W 3 (t - 0) 

tt 

X W 2 (t)w 2 (t) 

so that w 3 (t )is a functional of w 2 (t) w 2 (t). The operators should be evaluated 
in the asymptotic limit, which does not exist at this level because of particle 
streaming. This defect should be corrected by the inclusion of order fl effects 
in the solution for w 3 , but we may indicate the formal procedure . At the initial 
instant we have 


w 3 ( t r 0) 


3 

p 3 (t) w 3 (t = 0) ~ Y] £(i) w 2 (t 0)w 2 (t 0) 

tt ; t=0 


(26) 
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The left side denotes the measurable initial value for w 3 , from which the de- 
pendence onw 2 (t = 0) (also measurable) must be subtracted to yield the proper 
initial value for the term w 3 0 of equation 25. 

Due to the integrations of the operators 0(1, 3), o( 2, 3) on the initial value 
ofw 3 and the subsequent integration 0(1, 2) in the evaluation of the collision term 
for f, we may obtain an estimate of the effect of w 3 ( t = 0) on the evolution of f 
despite the unboundedness of the unintegrated quantities. 

The term w 3 (t - 0) contributes an order /3 2 term to the kinetic equation so 
that a rough estimate is satisfactory. By ignoring the non linear terms in 
Eq. (25), we may integrate to find this correction. 

(It) ( w 3 (t = 0)) = 0(1, 2) J drP 2 (t-r) 0(j , 3)P 3 (r)w 3 (t = 0) 

We do not evaluate this term here, for it has the general form of the col- 
lision integral produced by the initial value of g 2 of Section IV, and is formally 
order fi smaller. 

Henceforth we omit w 3 (t = 0) and use the convenient notation w 2 = (§fSf) c 
where the subscript denotes a continuous rather than a singular function. 


[ 1 2 3 /*t 

dtT + T( 1) + T(2) <Sf(X 1 , t) Sf(X 2 , t)) c = 0(j, 3) J o drP(\ i} t - r) 


0(i, t -r) P” 1 (X 4 , t-r)(Sf(-, t)Sf(X i , t)) c <Sf(-, t)8f(X 4 , t)> c (27) 
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Due to the sums 


2 3 

IX 

j=l i= 1 

and the two ways of distributing the remaining coordinates the right side con- 
tains 12 terms in (Sf8f) c (8f8f) c . 

We consider first the four terms containing both 0( 1, 3) and 0(2, 4). Upon 
evaluation we find that these terms have the form 

dkjdk' £(k, k') e ik '( r i‘ r 2) e ik '‘( r i‘ r 2) (8fSf | k) c <8f8f | k') c 

while all remaining terms have the general form 

dk e lk ( r i _r 2) J dk - £ ( k k 'j ( Sf Sf | k ) c (sfSf | k ') c 

The former terms correspond to an interaction between S f (X j) and Sf(X 2 ) over 
a range ) r x - r 2 | ~A. d , and will be neglected compared to the latter, which 
represent the interaction of Sf (X 2 ) and Sf(X 2 ) separately with a fluctuation 
background. This approximation, which is similar to the neglect of close en- 
counters in the equations for h n , causes a great simplification, for it follows 
that the behavior of {Xj } and {X 2 j separate as in Section V. It is convenient to 
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abbreviate notation by treating only the X x dependence, written | 8 f (Xj , t)). We 


have 


~ d -| f‘ r 

dT+T(l) I Sf(X x , t)> = 0(1,3) dr P(Xj, t) 0(1, 4) P" 1 (X 1 , r) P" 1 (X 4 , r) 

J Jo L 


{l Sf (Xj, t)> <Sf(X 3 , t) Sf(X 4 , t)> c + i Sf(X 4 , t)> <Sf(X„ t) Sf(X„ t)) c } 


+ P(X„ r) 0(3, 4) p- 1 (X 3 , r)p-‘(X 4 , t) {<Sf(X„ t) Sf(X 4 , t)> c | Sf(X„ t)> 

+ <Sf (x 3 , t) Sf(X 3 , t)) c | Sf(X 4 , t)U (28) 

The operators 0(1, 4) and 0(3, 4) create an electric field from the fluctuation 
density Sf (X 4 ). In order to proceed we assume that this field may be ex- 
pressed as 


f 


dk SE(k) e 1 


1 8 

Although this is consistent with the asymptotic behavior from linear theory, 
it does not necessarily give the proper connection to the initial state of the 
system, contrary to the procedures of Section III. Finally because the correlation 
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function (Sf (X.) Sf(X j )) is not a function of r . - r. except in the case of spatial 
homogeneity we use an adiabatic hypothesis to omit the dependence on slow 
spatial variation. 


it 


Sf(X t , ,)> 


0 ( 1 , 3 ) 




ik-r. 


q. JL 

m <?Vj 


x P- 1 (X,, T ) |(Sf (-k, v 3 , t) 8E(k, t)) c e- ik<r 3 


Sf(X 1 , t)>+ <Sf (X x , t) 8f(X 3 , t)) c 


8E(k, t) 


)} + P(X 3 - t) 


ik* r 3 d 
m dx. 


1 (X 3 , r) |(sf(-k, Vj, t) 8E(k, t)) c 


ik- r. 


x 


Sf(x 3 


t)> + <Sf(X 1 , t) Sf (x 3 , t)) I 8E(k, t) 



( 29 ) 


The P operator contains two terms, a flow operator and a shielding operator. 
Because evaluation of the latter requires once again an assumption about the 
behavior of the electric field we leave these terms for further study and ap- 
proximate 


P(X, r) e ik ‘ r P' 1 (X, r) by e~ Vr e ikr -37 e' 
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in Eq. (29). 


Ti * T t !>] I 0> 
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0(1, 3) 
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dk 




*T l( -r) yziSil 

dx , 


d dv i(T) alJ-ik-r, 


<9r j + dx 1 ' (?v,|j e ' 3 (Sf (-k, v 3 , t) 8E(k, t)) c | Sf(X 1 , t)) + (Sf(X lt t) 


Sf(X 3 , t)> c | SE(k, ‘» K 


ik-r,(-T) r r 3 ( r ) d * V 3( T ) d 


^ V 3 ^ r 3 + “'3 "’3 


dV~ ' a77h< Sf ( _k ’ v i* 0 


SE(k, t» c e' ik ' r3 | Sf(X 3 , t)) + <Sf(X 1 , t) Sf(X 3 , t)> | SE(k, t)V 


(30) 


The first term may be evaluated by carrying the r integration to infinity 



dk q 2 

(2tt) 3 m 2 


(SE(k, t) 8E(-k, 



ik" [r 1 (- T ) _r 1 l 
e L J 


(r) 

dx 1 


d dy i ( r ) 

dr t + dx t 



This term describes the diffusion of particle orbits due to background fluctua- 
tions in the electric field. This effect was first described by Dupree, 19 who based 
his theory on the Vlasov equation. 
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Before proceeding we discuss the consistency relations which modify the 
right side of Eq. (30). It should be noted that the desired consistency should be 
sought between the new equations for w n (Eq. (24)) where the behavior of each w n 
now differs from the /3 = 0 approximation because of interactions with w 2 . This 
modified behavior alters the P n which propagate the w n in time, which in turn 
leads to a correction of the "collision operator" relating dw n /dtto w m (t) w n _ m+2 (t). 

In the approximation that interactions between spaces X. , X. are neglected 
the consistency for the w n equations simplifies greatly because of the possibility 
of treating each "Sf" separately. The desired consistency among the equations 
for w n now reduces to self consistency of the | Sf) equation, as in the case of the 
Vlasov theory (Section B). 

We wish to modify P, which propagates | Sf)for /3 = 0, to absorb the correct- 
ions produced by the right side of Eq. (30). These corrections should alter in a 
self consistent way the collision operator relating d | Sf)/<?t to ($f(t) Sf(t))| Sf(t)). 
Since this appears very difficult to carry out we simply indicate how some of the 
consistency relations might be included. 

Thus if we modify the orbits and ignore the other order /3 corrections, then 
in Eq. (30) we replace (symbolically) the operator 


e -vr e ik-r 


dv 


,Vt 


by 


“V r ^ik*r 


(9 

d\ 


.V T 
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where 


V' 


d _d_ _d_ _d_ 

^ <?v ’ ’ dr d\ - ^v <9v 
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— -u n - q 

(9v 


+ D 


r <9r 
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dk 

(2tt) ; 


(SE(k, t) 8E(-k, 


Jo 


dr 


,-vV ;k-r 


_d_ 

d\ 


e v T e -ik ’ 


This correction has been considered in the limit of short time^ and long 
20 

time in other work. This orbit correction may be applied to any other 
resonance term to eliminate mathematical difficulties, but its effect may be 
small compared to the other order /3 terms in a particular physical problem. 

The fourth term may also be evaluated by carrying the r integration to 
infinity . 


q 2 d 

2 d v , 
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J (2T7) 3 J ( 2T7) 3 
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dv, 4rmq 


i(k + k ' ) 
Ik + k' I 2 


SE(k, t 



i o> k T 


x e 


ik-[r 3 (-T)-r 3 ] 


<?r 3 (t) d\ 3 (t) 

• ik' + 


d\ . 


d\. 


<9v. 


(§f (-k' , Vj, t) 8f [k' , v 3 , t)) 


This term represents a correction to the shielding of ] SE) due to non- 
uniformity in the background plasma. The principal effect is a correction to 
the dielectric function, leading to a change in the frequency co k . 
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The second and third terms describe effects similar to large angle scattering 


in a Boltzmann gas. The third 


dx 


r dk 

i J ( 2tt )' 


\M 


dv 3 47rnq 


i(k + k‘) 
(k + k') 2 


f 

» O 


dr e 


t ik-[r 3 (-T)-r 3 ] 


(Sf(-k, vj, t) SE(k, t)) c 


~dr 3 (r) 
<?v , 


ik ' + 


d hSll 

<3v. 


_d_ 
dx - 


8f(k', v 3 , t)) 


ik'-r. 


represents the correction to the electric field caused by the modified orbits of 
the fluctuation particles. The second cannot be evaluated as written in Eq. (30) 
because the time integral is not well defined at its upper limit. This may be 
corrected by including the diffusion of orbits, but the self consistent correction 
due to scattering is probably more important physically. In both terms the 
non-linearity must be taken into account in correcting the P operator, because 
d | 8 A/St depends on the velocity distribution of the correlation functions through 


i 


dv' 


dk' (8 f ( k ' 


v ) 8 f (-k ' , v ' )) 

' C 


We expect these scattering terms to be important for the description of acoustic 
phenomena where phase velocities are essentially constant over a range of wave 
numbers. 
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It is obvious that second order plasma theory is very complicated even 
when particle discreteness effects are omitted. It appears desirable to seek 
criteria for the relative importance of terms in a given physical situation, and 
means of approximating them, rather than a significantly better analytic treatment. 

A. Comparison with Quasilinear Theory 

In the analysis of quasilinear theory and mode coupling one derives an 
equation for a spatially averaged distribution f 0 (v, t) by solving the equations 
for the fluctuations (Fourier transformed) f h (v, t), E k (t). 

Thus the right side of the equation 

d q d f dk . . 

dt ~ m <9v ’ J (2n) 3 ^ V ’ (^) 

is to be found from the equations 


d q <? f Q q_ _d_ f dk' 

dt + ^ " v + m Mi dv + m dv " J (2n) 3 


0 (32) 


ik • E h 


4/t nq 



(33) 


The usual procedure is to solve Eq. (32) by perturbation theory, treating 
the convolution integral as small. Thus the first approximation arises from 
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the solution of the linearized Vlasov equation, while the corrections may be 
written as the sum of two series: 1) Terms proportional to df 0 /dv which we 
call shielding terms and shall discuss later. 2) Terms containing 
e -ik- yt f t = 0) which are usually called the "initial value terms." We shall 
discuss these terms first. Dropping the shielding term from Eq. (32), and per- 
forming a power series expansion in E, we find 

fk.i.v. = e ' lk ' v ' M*. t = 0) 


p 1 

k, i. v 
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dr 


a -ik *v (t-r ) 
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( 9 v 
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dk* 

(2 rry 


E k' (t - 0 ) 


f k-k J ,i.v. ( V ’ T ) 


where the dependence on f(t - 0) arises from the convolution over f k-k ' • 
Because 



^k + m dv 


i 


dk' 

(2 V ) 3 


"k' 


^k-k' 


is the transform of 


(dt + v * V + m SE ‘ dv) 8f 

which may be identified with particle orbits, we draw the following conclusion: 
Neglecting initial value terms is equivalent to neglecting the effect of electric 
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fields on the orbits (or statistical orbits) of particles in the fluctuations, and 
the modification of the electric fields due to these particle deflections. 

As an example we consider f k 2 i v from the perturbation analysis and 
average over the electric fields in the random phase approximation. 


f 2 

r k,i.v. 


Jo 


dr 


a -ik -v(t-r ) 


J_ f dk' 

(9v J (2tt) 3 
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ik* Vr' 


f k (v, t - 0) 


The quantity is equal to the first correction (in D r ,D v ) to the solution of 


r<? d e d „ 

|dt + v * V " dv ’ D v * dv “ dv ' D r * drj Sf “ 0 

where D r and D v are given without the consistency requirement, and with upper 
limit t in the time integration. By dropping initial value terms we omit this 
correction. 

We next consider the effect of the perturbation treatment on the shielding 
terms. A simple way to estimate some of these effects is to expand (inD r , D v ) 
the equation 


fd \ d d d 

^dt + ik ' v j fk ~ dv * ^v dv ~~ dv 


ik f k + m E k 


df, 
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0 
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The lowest order term 


fk°, = fh-il-')" (e- 1 “‘ , -e-*' vt ) E„ 


(t = 0) 


vanes as 



E k 


df Q 

dv 


in the resonance region co k = kv. We may estimate the fastest growing correc- 
tion from the term in D v , and overestimate the time of validity of the expansion 
by treating D v as constant. 
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By comparing f k ' s and f k ° s we see that the perturbation series breaks down 21 
in the resonance region for t ^ (k 2 D v )“ 1/3 . From the equation 


^£0 i_ 

fit _ dv ’ ^v ' dv 
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we may estimate the time for f 0 to change as t ~ v 2 /D y . Thus the perturbation 
expansion may break down before f 0 changes significantly. 

Although this breakdown need not affect the electric fields strongly (the 
secularities vanish on velocity integration) it should modify the behavior of f 0 
significantly. We believe that corrections to the f 0 equation (through f k ) should 
be included if the mode coupling analysis is to be considered valid. 

B. Comparison with Dupree's Theory 

In Dupree's perturbation theory for plasma turbulence the starting point 
is the Vlasov equation together with Poisson's Equation (33). 

/ d „ q d\ 

\ft + v ‘ ^ + m E ' Tv) f ( r ’ v - = 0 (34) 

The electric field and distribution function are then expanded in Fourier 
series in space (k) followed by an additional Fourier series in the phases 
( At ’ At ' * * * ) relative to the arbitrary initial phases of the electric field. 

E(r, t) = Efc (O ex p(ik • r + i/S k ) 

k 


00 
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Before proceeding we point out that Dupree's statement that at the initial 


instant 


F ( n k , n k ' , • • •) - 0 unless n k , n k < , n k " • • • = 0 

is not consistent with Poisson's equation, so that the electric field is not uniquely 

determined. The statement regarding f(t = 0) is essentially that of dropping the 

initial value terms of quasilinear theory,^ except that the diffusion correction 

D v to the orbits is produced by the statistical fields. The consequences of this 

19 

approach have been discussed by Orszog and Kraichnan, who have noted the 
similarity to the stochastic acceleration problem. 

Instead of discussing the difficulties of Dupree's approach, we give a simple 
procedure based on the Vlasov equation which substantially duplicates the theory 
based on the equations for w , while permitting a comparison with Dupree's 
result. 

We may ensemble average Eq. (34) in any convenient way to find 

{it + v • V )( f > + rn ( E > • ~k < f > = "miv'( SfSE ) < 35 > 

The difference quantity f - (f) = Sf then satisfies the equation 

(£ + * • 7 ♦ 1< E > • 77 ) «f = ~f>m 77- [ 8fsE - <S f SE>] <3G) 
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where we have inserted /3 in order to follow the procedure described in Sections II 
and HI. We wish to obtain a statistical estimate of the order fi terms by using a 
self consistent perturbation theory. From a power series expansion in /3 we find 


Sf° = P(t) Sf(t = 0) ; 


8E° 


P(E, t) 8f(t = 0) 



SE(k, t = 0) 


8 f 


“ m J Q drP(t-T)^- • [Sf(r> SE(t) - <8 f (r) SE(r)>] 


8E 1 



Sf 1 (r' t v\ t) 


By using the /3 = 0 solutions we may write 8 f 1 ( t ) and SE 1 (t ) in terms of 8 f 0 ( t ) 
and SE° ( t ). Then multiplying 8 f 1 by SE and SE 1 by 8 f and dropping sub- 
scripts yields 
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We now average in all ways so as to obtain pair correlations, and add to obtain 


all statistical contributions to the right side of Eq. (36). 
statistical^ _ ^ d f‘ f dk , 

collisions / m 2 ^ J 0 J ( 2n ) 3 6 


P(X, T)e ik - r ^ r P" 1 (X, t) 



x {<SE(k, t) SE(-k, t)) e~ ik r 5f(X, t) + <Sf(X, t) SE(r, t>) SE(k, t)j 

-TT J*' |7T7q P<X',T)e‘K- r '^ T p-' ( X',T) 

x j(SE(k, t) Sf (— k , v, t>) e _ ik ' r S f ( X ' , t) + (§f(X' , t) 8f(X, t)> SE(k, t ) j 

Although we are unable to motivate the above procedures, the result is identical 
to Eq. (29). As discussed previously one of these terms (modified by the con- 
sistency relation) is the diffusion correction to the orbits first obtained by 
Dupree. 

VI. Boltzmann Gas 

In this section we use the BBGKY equations for g n (X, • • • X n ), the n par- 
ticle correlation function - f 'j . The force field F is produced externally 
and will be considered independent of time. Although some results may be 
obtained more directly by considering the n particle distribution f n , we follow 
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I 


the procedure described previously. g n satisfies the equation 



4> is the potential, and fi is the ordering parameter. The sum 

m= 1 xt j 

includes all possible ways of distributing the n coordinates such that each ap- 
pears once, and the i, j sum runs over all values which take one coordinate 
index from each correlation function. 

Instead of attempting to write the full solution for /3 = 0, we consider the 
effect of the source terms for g n as successive corrections, and use a second 
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subscript to indicate the order of approximation of the source terms. Thus g n 0 
(no source term) is equal to e -Hnt g n (Xj, • • • X n , t = o) , where e~ Hnt runs the 
particles backward on their orbits for time t. By using the g m 0 solutions in 
the source terms we may find g n j . Since 


r. V i ' <?V. X i’ •") 8n-n.(‘“- X j- ’“l = [ H n( X l- ' * ’ ’ X n) 


i . J 


^ H„(X„ X b ,X t , ■•••)- H„..(X p , X,, X, '••)]g.(X„X b , X c . •••) g„.„(X p , X q ,X r ,--) 


we may write (with n > 2 , j3 - 0) . 



= 0 



= 0 


s n , 1 


e 


n- l 


-H t 


- 1 n ~ t 

g m ,o (t = °) g n - m ,o (t = °) " 2Z gm '° (t) gn ‘ m '° (t) + e " g " (t = 0) 

= 1 m ~ 1 

(38) 


Note that the sum on m contains all distinct ways of distributing the n coordinates 
between g m and g n _ m . We use the relation g m 0 (t = 0) = e 11 ™ 1 g m 0 (t) and then 
treat the g m 0 as the exact functions, thereby throwing the approximation over 
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to the operator . 


*„.!<») = £^{e" H "'e""' e""-"'- l} g„(t)g n _„(t 


) + e Hnt 8 n (* " °) 


The limit t-c° on the operator does not exist unless boundary conditions at 
r = oo are specified on the functions g m (including gj - f !), for this operator 
brings in particles from greater distances as t increases. This should not 
occur because the orbits are altered by statistical encounters with particles in 
the space (r n+ j , v n+ ^ . However these effects are excluded by the ordering of 
terms in the hierarchy, so that order /3 effects should be included to give a 
proper result. Nevertheless g 2t i (/? = 0) may be used to obtain the first ap- 
proximation to the collision integral. 

A. Boltzmann Equation 

For n = 2 we have 


B 2 


-H,t 


g 2 (t - °) + 


{' 


■H 0 t H,t H , t 


*) f ( x 2 . ») 


(39) 


We substitute this result into the equation for f , and discard the initial value 
term g 2 (t= 0) , which will vanish for long times if order fi effects are included 
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in the g 2 equation. 


d F d 

Jt + V 1 • v i + m • <9 Vl 


n Jdr 2 dv 2 Vj <f >1 2 • dv ^ [e 2 e 1 e 1 l] f (X 1( t) 


f (X 2 , t) + n J*dr 2 dv 2 Vj 2 • ^ f (Xj , t) f (X 2 , t) + n [1 J*dr 2 dv - 


2^1 ^12 


• ^7 [g2 - g 2 (^ - °)] ( 40 > 

The term g 2 ' g 2 C/ 3 = °) is formally order /3, so the corresponding term in the 
equation is order /3 2 and may be neglected. We cut off t in the operator at some 
large time T, use the consequent relation 


H 2 g 2 = -v ,< P 


12 


(*-£) 


ff 


and the fact that 


JdX 2 V 1 0 12 • ^ - JdX 2 [H 2 -H, (1)-H 1 (2)] 


to write 


,1 F 

v l ,V i + m - <?v 


-If - nJdr I dv I [v,.V, t v 2 .V 2 ][/ B2T e , ' T e H,T -l]f(’‘,.t)fWt) 
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Although Bogoluibov-*- obtained essentially this result he evaluated it in- 
correctly. The result is improper as written because the correlation extend 
to large distances |vj - v 2 | T, or q° in Bogoluibov's derivation. Since interactions 
with a statistically distributed third particle (an order /3 effect) will destroy the 
correlation we cut off the integral at a distance k somewhat less than a mean free 
path. At distances of this order the collisional correction to g 2 , and to f (since 

-H.t 

f does not evolve as e f(t - O)for a mean free time) must be taken into ac- 
count by the second order theory. 

We first change the variable of integration tor = r 2 -r, = (b, 4>, z) where 
z is parallel to v 2 - v j, and use the fact that 



J 


dXg 


0 


n 



bdbdc^dz 



(9 

d z 



HjT HjT 

e e 



X f (r,, Vj, t) f (rj + r, v 2 , t) 


The external force affects only the motion of the center of mass and drops 
out of the operators in the collision integral. The z integral is cut off as 
specified above. 




J 


bdbd (f) | v 






y 


*) 


- f ( r i- v i) f ( r i + X, v 2 )} 
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Here r^Vj', r 2 ', v 2 scatter in the absence of an external field to produce r 1# 
Vj, rj + A., v 2 . We expand r and r 2 ' above r ,, keeping terms in Vf , and find 
by evaluation 


Ar 


l 


^( v i' ~ v i) 

I v l - V 2 I 


so that the resultant equation is 



d f F 

dt + v i • V i f + m 


• ~ n jdv 2 Jbdbd<£ \y 1 ~ v 2 | [f(r,, v/) f(r lt v 2 ') 


(r,, v,) f (r j , v 2 )] + n\ f dv 2 j" bdbd p 


<?f (r P v i') 

( v i' - v i) • f ( r i- v 2) 


+ (v 2 ‘ - V,) • f (r j , v,') 


df ( r i> v 2 ') 

dr 1 


- n K 



bdbd0 (v 


- v 


2) 


' f ( r i- v i) 


<5 f (r 1 , v 2 ) 
( 41 ) 


The first term on the right is the usual Boltzmann integral, and the cor- 
rections due to a spatial gradient may be identified with the orbits of the scat- 
tering particles. The first represents a correction to the number scattered 
into r j, vj because the number of collisions at a distance k is not generally the 


50 



same as the number at rj. Since particles scattering into r I( Vj suffered col- 


lisions at an earlier time, for consistency the number scattered out of r ls v, 
should be evaluated at this earlier time. This gives rise to the second term 
in Vf, a correction to the number scattered out of r 1( v 1( because the target 
particle density must be evaluated at r j + \ . 

The collision integral conserves local number density, but it contributes 
to a momentum and energy flux. 

j” dv j x(v j) [collision integral] = nA.Vj J dv j v j y (v x ) J dv 2 jf (r j , Vj) 

x f(r„ v 2 ) - f(v lt v,') f (r lt v 2 ')j 

where 


x(vj) - jnmv j , 

We have not pointed out several minor difficulties in the derivation of the 
kinetic equation, as their resolution should be based on the inclusion of three 
particle effects. 


i- X 

2 nmv ! J 
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B. Higher Order Effects 


As in Part A we substitute each g n+ i i into the order /3 term of the equation 
for g n and calculate the eollisional correction to behavior of g n . 





tT 


[g n+ i.i-g n +i] (42) 


The collision term is evaluated as before. 


dt 


”1 n_ ^ r 

H „l *.-£]2Z' 7 ‘ **j - ^7 •;.(••• x i-) «»-. ( - • x i • ■ " J dX »* ■ 

m- 1 i , j 


n+ 1 “1 

1=1 



H T H 

e m e n 


n + 1 1 



- -p/3 



^+i 


— — r 
dv. I^n+l 

i L 



Each integral over d/dr n + x should be cut off consistent with the collision oper- 
ator in the next higher equation. We are unable to perform this operation at 
present, but progress is still possible for n = 2 0 
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For n — 2 we may write 



§2 


■ 3?7 f f -s k^lv’J £ 

j^j= i " i=l permutations 


X L-V e H 2 T e"‘ T -l}g 2 (X.,X.,t) f(X„.t) = -n/3 J* dX 3 [h 3 - H 2 (X, , X 2 ] “H,(X 3 )] 

x |g 3 2 ~ §3,i} + n /3 2 j dX 3 [H, - H 2 (X,, X 2 ) - Hj (X 3 )J |g 3 .2 _ g 3 } (43) 

where g 3 2 includes the second approximation to the source terms of the g 3 
equation. The difference between g 3 and g 3 2 is regarded as order /3 and will 
be neglected here. (In general g n n _ t is a functional of f , and the difference 
g n - g n _ j should be treated as one order smaller, i.e. g n - g nn _i = 0(/?g n ))- We 
may calculate g 3 2 to order 1 by using the first approximation g 2>1 in the source 
term. 



§ 3,2 


♦E L ^ 

i , j permutations 


d 

d\ . §2, 1 f 


y [h 3 -h 2 -hj g21 f 


53 




In calculating the difference £(g 3 2 ~ g 3> i) to order /3 we may approximate 
83.1 (62) to order one by g 3>1 (g 2 x ) 


x f(X i , t) f(X.. f) f(X k . t) 


Taking the difference we find 


- 83 ,, = £ e“ 3 ‘ e H2<X > ' X j )l e »,<X k )* _ *-.• .«.(*,)• 


« e H,<Xi) ' e“‘ (Xk) ' + t) f(X r t) f(\, t) 


We substitute this result into Eq. ( 43 ) to find the equation for g 2 , correct 
to order /?. 

[W + ** 2 ] g 2 - V i j ’ f f - n f dX 3 J^[ v i ' v i] 2^ | e HsT eH2T e ” lT 

i + j J i= 1 P t 


X g 2 f - 


- S J d X 3 g [vi .V i] /^e-Ve^e"' 


T -H,T H,T H.T H,T 

- 2e 3 e 1 e 1 e 1 + 1> fff 


( 44 ) 
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The integrals should be cut off by considering the effect of a statistically dis- 
tributed fourth particle. By initial hypothesis g 2 should depend on this cutoff 
in a weak way (order / 3 2 ). As before g 2 should be obtained in terms of an 
arbitrary f, after which the dependence on f at earlier times should be elimin- 
ated using the formal solution for the Boltzmann equation (Eq. (41)). The latter 
requirement may be avoided if we consider a system sufficiently near equilibrium 
so that the change of f due to flow is canceled to order (i by the effects of col- 
lisions. Even in this (transport) regime the solution of the g 2 equation appears 
difficult, for the collisional shielding terms cannot be found by expansion. 

VII. Plasma Equilibrium 

For the study of thermal equilibrium the cluster expansion discussed 
recently by Ramanathan 23 is most suitable. In equilibrium the one particle 
distribution is simply f (X x ) = f m (v 1 ), the Maxwellian velocity distribution. 

The pair and triplet correlations a. . (r. - , a. jk fr { , r . , r k ^ are defined 

by 

f 2 = f m( v l) f m ( V 2 ) l 1 +a 12) 

f 3 = f n( V l) f m( V 2) f m ( V 3 ) ( 1 + a l 2 ) ( 1+a l 3 ) ( 1 + a 2 3 ) ( 1+a 12 3 ) 

The pair correlation is related to that of the BBGKY by g 2 ~ f m (v j) f m ( v 2 ) a 1 2 
while the triplet correlations are quite different. For either hierarchy of equa- 
tions the usual expansion procedures are not uniformly valid; they break down 
for large r. .. In this section we find a convergent order 0 correction to a i 2 - 
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The function a 12 satisfies the equation 


C( 1 ) a j 2 


- V, 0 


12 KT 



+ a 



+ 0(/3 2 ) + 


where the operator C is defined by 


C(i)0(X.,X., ...) 


KT J 


V. 0 + kt dr j< v. 4 >. ti , 0(r i 


and cpis the potential 4 > i . = -q 2 /|r 4 - r ^ |. The triplet correlation a ijk 


the equation 


c (i) a ijk 


.if 

KT J 


dr 4 V. 4>. 4 a j4 a k4 + 0(/3) + 


and all /3 = 0 equations for n > 2 have the form 


C ( 1 ) a i j • • • n - “ KT j dr n+l V i ^i.n+1 ^ 


no dependence'! 


on r . 


a 12 a 23 


• • • (45) 


) 


satisfies 


(46) 
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We now define which is the zero order approximation to the pair correlation, 

by 

c (i) Xij = " W V i ^ij Xij (°°) = 0 

In a plasma 




KT r. - r 

' 1 J 




where A. , 2 = KT/ 477 nq 2 We see at once that Eq. (46) has the solution 


a. .. 

i jk 


n f 


dr 4 x i4 a j4 a k4 


We substitute this result into Eq. (45), with i = 3, and use the approximation 


'23 


X 23 bi the term a 12 a 23 , to find the equation for a 12 correct to order fi. 


C( 1) a 12 


KT v J 0 1 2 + fi 


V 1 Xl2 a 12 + n 1 

m 

dr 3 Vj X 13 a i3 a 2 3 

J 

- 


(47) 


If we use perturbation theory we substitute a. t . - ■ on the right and use the 

integral 


C (i) a i... = •••) 4 a i... 


^( r i- ’••) 


n | dr i' X ifi ' 0( r i'* ‘ - 
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then we find at once 


1 

a 12 “ X 1 2 + 2 I ( 1 ) 

24 

which is the result obtained previously by Rostoker and O'Neil for r 12 ^X d . 
This result breaks down for r l2 ~q 2 /KT, where a 12 = e* 12 - 1, and for r 12 >> \ d , 
In the latter case we approximate a t 2 and a j 3 by x in the order fi terms of 
Eq. (47), to find 


X,V 




dr 3 Xm X23 


C( 1 ) a i 2 



V 1 Xl 2 3 a 


23 


1 

- KT Vj 0 12 



(48) 


We use a Fourier transform to solve Eq. (48). 


-12 


(k) 


n 



i i ( q 2 

) i . 

A ~ ik ^d/ 2 \ 



k 2 \ d 2 " 4 

' kX d 2 ln 

V 1 + ik ^d/ 2 / 


1 

. 1 if' 

q 2 \ 1 

A - ik ^ d / 2 \ 

k 2 x d 2 " 4 V 

KT/ k X d 2 

ln \1 + ik X d / 2 / 


The inverse transformation is evaluated by analytic continuation. The result 
comes from the pole at 


k 


o 



8 KT \ d 2 


In 3 
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where the contribution from a branch cut starting at k = 2ik d may be neglected 


for large r , 2 . 


' 1 2 


(r) = 


- k 


KTr 


o r 


8 KT \ 


In 3 


( 49 ) 


6KT\ 


It is apparent that the asymptotic behavior of a , 2 remains undetermined. 


VIII. Conclusion 

Although the procedures described in this paper appear adequate for the 
calculation of a reasonably accurate (first order) collision term, it may be de- 
sirable to recast the theory in less mathematical form for the calculation of 
higher order corrections. Since the solution of the first order theory is necessary 
for the calculation of higher order terms, we must consider the development of 
procedures for solving kinetic equations as a primary goal. For the important 
case of asymptotic solutions (in the domain of transport theory) the knowledge 
of first order solutions may not be necessary, but it is not clear that second 
order corrections will be important in this regime. 
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